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^ ' Abstract 

^^ I We analyze the angular dependence of the irreversible magneti- 

^^ ■ zation of YBa2Cu307 crystals with columnar defects inclined from 

^-^ , the c-axis. At high fields a sharp maximum centered at the tracks' 

Q^ I direction is observed. At low fields we identify a lock-in phase char- 

^^ ' acterized by an angle-independent pinning strength and observe an 

"t^ ■ angular shift of the peak towards the c-axis that originates in the ma- 

terial anisotropy. The interplay among columnar defects, twins and 
ab-planes generates a variety of staircase structures. We show that 
correlated pinning dominates for all field orientations. 



Q ! A difficult aspect of the study of vortex dynamics in HTSC in the presence 

of correlated disorder is the determination of flux structures for applied flelds 
'kp I tilted with respect to the pinning potential. As 3D vortex conflgurations 

;_( ■ cannot be directly observed, our knowledge is mostly based on the analysis of 

the angular dependence of magnetization, susceptibility or transport data[|l], 
I, I, I, I, I, 0,1, 1,0. 



According to theoretical modelslTJ, 0, when the angle between the ap 



plied fleld H and the defects is smaller than the lock-in angle (pi vortices 

1 



remain locked into the defects thus producing a transverse Meissner effect. 
For tilt angles larger than ipi and smaller than a trapping angle ifT-, vortices 
form staircases with segments pinned into different defects and connected 
by unpinned or weakly pinned kinks. Beyond ifx-, vortices will be straight 
and take the direction of the applied field, thus being unaffected by the cor- 
related nature of the pinning. In principle, this picture should apply with 
minor differences to twins, columnar defects and intrinsic pinning jT^. 

Many experiments have confirmed the directional pinning due to colum- 
nar defects, twins and Cu-0 planesfll, ||, |], H, ||, P, |^, P, ||, |T^. Evidence 



for a locked-in phase arises from the observation of the transverse Meissner 
effect 0, but a quantitative determination of ipL{H,T) for columnar defects 
had not been done until now. The introduction of columnar defects inclined 
with respect to the c-axis has been used[|l], |^, ^ ||, [1^] to discriminate their 
pinning effects from those due to twin boundaries, and from anisotropy ef- 
fects. However, the combined effect of the various correlated structures on 
the vortex configurations remain largely unexplored. 

In this work we report studies of the vortex pinning in YBa2Cu307 crys- 
tals with inclined columnar defects, for the whole range of field orientations. 
This allows us to determine the misalignment between the applied and in- 
ternal fields due to anisotropy, as well as to identify the angular range of 
influence of each correlated pinning structure. We present the first determi- 
nation of the lock-in angle of tracks using irreversible magnetization. 

The crystal used in this study was grown by the self flux method [|T^] , 
and has dimensions ~ 200 x 600 x 8.5fim^. Columnar defects at an angle 
0£) ~ 32° from the c-axis and a density corresponding to a matching field 
Bq, = 3T were introduced by irradiation with 315 MeV Au'^^^ ions at the 
Tandar accelerator (Buenos Aires, Argentina). 

DC magnetization M was measured in a Quantum Design SQUID mag- 
netometer with two sets of pick up coils, and both the longitudinal {Mi, 
parallel to H) and transverse (M^, perpendicular to H) components were 
recorded. The sample could be rotated in situ around an axis perpendicular 
to H using a home-made device. The angle B between the normal to the 
crystal (that coincides with the c-axis) and H was determined with absolute 
accuracy ~ 1°, and relative variations between adjacent angles better than 
0.2°. The details of the procedure are described elsewhere |I1 . 



Magnetization loops Mi{H) and Mt{H) were recorded at fixed T and 0. 
Sample was then rotated, warmed up above Tc and cooled down in zero field 



to start a new run. We use the hysteresis widths AM^ (H) and AMj (H) to 



calculate the modulus Mj = ^J ^Mf + AM^^ and direction of the irreversible 
magnetization vector Mj. It is known that in thin samples Mj is normal to 
the surface due to geometrical constrains |P^, IT^l, except above a critical 



angle (of ~ 87° for our crystal). We confirmed that Mj || c within 1°, for all 
e < 85°. 

From now on we analyze the modulus Mi as a function of T, H and B. 
Figure 1 shows Mj(B) at 60 and 7QK and several values of H . According 
to the Bean model, Mj oc J, where the screening current density J is lower 
than Jc due to thermal relaxation. The geometrical factor Mj/ J depends on 
0, but is almost constant for B below the critical angle. 

The most obvious feature of Fig. 1 is the asymmetry with respect to 
the c-axis, which is due to the uniaxial pinning of the inclined tracks. At 
high fields [H > IT) we observe a large peak in the direction of the tracks 
Od ~ 32°. For H < IT the peak becomes broader and progressively shifts 
away from the tracks in the direction of the c-axis as H decreases. The shift 
decreases with increasing T as shown in figure 2, where the angle 0max of 
the maximum in Mj is plotted as a function of H for three temperatures. 
The inset of figure 2 shows a blow-up of the data of fig. 1 ioi H = OAT and 
T = 60K. This curve exhibits the second main characteristic of the low field 
results, namely the existence of a plateau in Mj (6) (We define 0max as the 
center of the plateau). 

We first discuss the origin of the shift. Maximum pinning is expected 
to occur when the tracks are aligned with the direction that the vortices 
would have in the absence of pinning. For an anisotropic material, such 
direction does not coincide with H. If Gb is the angle between the equilib- 
rium induction field B (which represents the vortex direction) and the c-axis, 
minimization of the free energy for H^^ -C if -C H^2 gives ||12|| 



sm (Bb — B) ?^ — — — -— — In 



2H\nK £(Bb) \ B ) ' ' 

where Hc2{Qb) = -^c2/^(©b)- Here H^^ and H^2 are the lower and upper c- 

axis critical fields, e is the anisotropy and e (9) = f cos^ 6 + e'^ sin^ ^ j . For 
e < 1 vortices tilt towards the ab plane. When B = Qd we have B^ > Bd 
and the optimum pinning situation is not satisfied. Instead, maximum Mi 
occurs at the vortex-track alignment condition Bb = B^. This corresponds 



to an applied field angle Omax < ©d that can be calculated from Eq. [I| by 
setting Qb = ©d ~ 32°. (In this picture the peak cannot occur at 9 < 0, 
thus the negative values of 6max given by Eq. [I| at low H are unphysical, 
and Omax ^ as H —^ 0). 

The solid lines in fig. 2 are fits to Eq. |l]with fixed parameters [p!2| e = 1/7 



and H^2{T) = 1.6T/K x (T^ — T) (the fits are not very sensitive to any of 
them). Using H^^,{T)/2\nK = ^/SttXHT) and Xl{T) ^ A2,(0)(1-T/T,)-i, 
we obtain a good fit to the data as a function of field and temperature by 
setting only one free parameter, Xabi^) ~ 500 A. Although this value is 
significantly smaller than the accepted value [Q (~ 1400 A), we neverthe- 
less consider that this simple model captures the basic physics. We note 
that Zhukov et al. ||^ have reported lock-in angles for twin boundaries in 
YBa2Cu307 crystals that imply an H^^ about 5 times larger than the usual 
values, a result suggestively similar to our case. 

We now return to the plateau seen in the inset of figure 2. The constancy 
of Mi (6) indicates that the pinning energy remains constant and equal to 
the value at the alignment condition B^ = Gd- This behavior is a fingerprint 



of the lock-in phase |Tl|. The extension of the plateau in the H — Q plane at 
QOK (determined with accuracy ~ 1°) is shown as bars in Fig. 2. Its width 
decreases approximately as H~^, as expected |1TT|] for ipi, and for H > IT it 



becomes undetectable with our resolution. The decrease of Mi at the edges of 
the plateau is sharp, a result consistent with the appearance of kinks, which 
not only reduce Jc but also produce a faster relaxation. 

When |0^ — 0£)| > (fL vortices form staircases. Two questions arise here. 
First, which is the direction of the kinks that connect the pinned portions of 
the vortices? Second, do we observe evidence for a trapping angle (^5^? 

For > 0max, there is a wide angular range in Fig. 1 in which Mj (-1-0) > 
Mi (—0) for all H, i.e., pinning is stronger when H is closer to the tracks than 
in the crystallographically equivalent configuration in the opposite side. This 
asymmetry demonstrates that at the angle -f vortices form staircases, with 
segments trapped in the tracks. For < 0max we again observe asymmetry. 
Mi (0) crosses = with positive slope, indicating that pinning decreases 
as H is tilted away from the tracks. We can conclude that staircases extend 
at least beyond the c-axis into the < region. 

The angle 6k between the kinks and the c-axis can be calculated by mini- 
mization of the free energy [|12|. For simplicity, we consider the case H ^ H^^, 



where 0^ = and the problem reduces to calculate the energy of one vortex. 



as the other terms in the free energy are the same for all configurations]^, |1^ . 
If Lp is the length of a pinned segment, and Lk the length of the kink 
(see sketch in figure 4), the energy is -E oc Lptp^Qn) + Lkef{9k), where 
e/ (^fe) ~ ^0^ iPk) [liiK + 0.5] and Cp (©d) ~ £^o^ (©-d) [1^^^ + (^t\ are the line 
energy for free and pinned vortices respectively, Eq is the vortex energy scale 
and at < 0.5 parametrizes the core pinning energy due to the tracks (smaller 
at implies stronger pinning). Minimizing E with respect to 9^ we obtain two 
orientations, 6^^ for 6 < 6d and 0~l for 6 > Qd- 

As the tracks are inclined, 9^ and 
gles are independent of 6. As |0 — ©^ 
Lp decreases and the number of kinks increases, consequently the pinning 
energy lowers. This accounts for an Mi that decreases as we move away from 
the tracks. In particular, for 9 = ^^ vortices become straight {Lp = 0), thus 



et 



are different. However, those an- 



increases, 9^ remain constant while 



^T 



9i-Q 



k -"^D 



are the trapping angles in both directions. In general 9f^ 
must be obtained numerically, but for etan^^. ^ 1 and etan©^) ^ 1 we 
obtain 




tan9*«tane.±^,/j;^^^^^ (2) 

Eq. 1^ adequately describes the main features of the asymmetric region in 
Fig. 1, and for Qd = it coincides with the usual estimates [ITT], O of v?r- 



There is, however, an important missing ingredient in the standard de- 
scription presented above, namely the existence of twins and Cu-0 layers, 
which are additional sources of correlated pinning. This raises the possibility 
that vortices may simultaneously adjust to more than one of them, forming 
different types of staircases. 

Pinning by twin boundaries is visible in figure 1 as an additional peak 
centered at the c-axis for H = 2T and T = 60K. A blow-up of that peak 
is shown in the inset. We observe this maximum for H > IT. The width 
of this peak, ~ 5°, is in the typical range of reported trapping angles for 
twins 0, |, H, |, ||. On the other hand, the fact that the peak is mounted 
over an inclined background implies that vortices are also trapped by the 
tracks. Thus, vortices in this angular range contain segments both in the 
tracks and in the twins. These two types of segments are enough to build up 
the staircases for > 0, but for B < a third group of inclined kinks with 
9k < must exist in order to have vortices parallel to H. 

Another fact to be considered is that there is an angle Qsym (which is only 



weakly dependent on H) beyond which Mj (B) recovers the symmetry with 
respect to the c-axis. This is illustrated in Fig. 3, where Mj data for — |B| 
was reflected along the c-axis and superimposed to the results for + |6|. 

One possible interpretation is that for > Qsym staircases disappear, 
i.e., that 0^^^ = ^t ^^^ ^^ ^^^ determining ip^ = Qgym — ©d- However, this 
is inconsistent with our experimental results. Indeed, (p^ should decrease 
with T, and this decrease should be particularly strong above the depinning 



temperature |T6| Tdp ~ 4:0K due to the reduction of the pinning energy by 
entropic effects [|lll. This expectation is in sharp contrast with the observed 
increase of Qsym with T, which is shown in Figure 4 for H = 2T. Thus, 
the interpretation of Qsym as a measure of the trapping angle is ruled out. 
Moreover, if in a certain angular range vortices were not forming staircases, 
pinning could be described by a scalar disorder, then at high fields Mj (0) 
should follow the anisotropy scaling law[]T7[ Mj {H, 0) = Mj {e (0) H). Con- 
sistently, we do not observe such scaling in any angular range. 

Our alternative interpretation is that, at large 0, the kinks become 
trapped by the a6-planes. This idea has been used by Hardy et al. [0 
to explain that the Jc at low T in the very anisotropic Bi and Tl compounds 
with tracks at 0d = 45° was the same for either = 45° or = —45°. Our 
situation is different, as we are comparing two kinked configurations. 

We first note that, according to Eq. 0, 5^ cannot be exactly 90° for finite 
e, thus the intrinsic pinning must be incorporated into the model by assign- 
ing a lower energy to kinks in the a5-planes. Vortices may now form struc- 
tures consisting of segments trapped in the columns connected by segments 
trapped in the a6-planes, or alternatively an inclined kink may transform into 
a staircase of smaller kinks connecting segments in the planes (see sketches in 
figure 4). We should now compare the energy of the new configurations with 
that containing kinks at angles O/^ . This is equivalent to figure out whether 
the kinks at 6'^ lay within the trapping regime for the planes or not. The 
problem with this analysis is that, as 9^ are independent of 0, one of the 
two possibilities (either inclined or trapped kinks), will be the most favor- 
able for all 0. Thus, this picture alone cannot explain the crossover from an 
asymmetric to a symmetric regime in Mi (0). 

The key additional concept in this scenario is the dispersion in the pin- 
ning energy. The angles 9^ depend on the pinning strength of the adjacent 
tracks {at in Eq. |), thus dispersion in at implies dispersion in 9^. As 
increases, it becomes larger than the smaller 6'^ 's (that connect the weaker 



defects) and the corresponding kinks disappear. The vortices involved, how- 
ever, do not become straight, but remain trapped by stronger pins connected 
by longer kinks with larger 9^. This process goes on as G grows: the weaker 
tracks progressively become unnefective as the "local" 9^ is exceeded, and the 
distribution of 9^ shifts towards the afe-planes. When a particular kink falls 
within the trapping angle of the planes, a switch to the pinned-kink structure 
occurs. In this picture, the gradual crossover to the symmetric regime takes 
place when most of the remaining kinks are pinned by the planes. 

If kinks become locked, the total length of a vortex that is trapped in- 
side columnar defects is the total length of a track, independent of 6, and 
the total length of the kinks is ex tan(|B ± Od|) for field orientations ±G 
respectively. As |0| grows, the relative difference between the line energy 
in both orientations decreases, an effect that is reinforced by the small line 
energy of the kinks in the a6-planes. If kinks are not locked but rather form 
staircases, taking into account that the trapping angle for the ab-planes is 
smallQ] (~ 5°), the same argument still applies to a good approximation. 
The temperature dependence of Qsym is now easily explained by a faster de- 
crease of the pinning of the a6-planes with T as compared to the columnar 
defects. Additional evidence in support of our description comes from recent 
transport measurements in twinned YBa2Cu307 crystals, which show that 
in the liquid phase vortices remain correlated along the c-axis for all field 
orientations^^^, suggesting that they are composed solely of segments in the 



twins and in the a6-planes. 

In summary, we have shown that the combined effect of the three sources 
of correlated pinning must be taken into account to describe the vortex struc- 
ture in samples with inclined columnar defects. We demonstrate that the 
lock-in phase exhibits an angle independent pinning strength, and show the 
decrease of the lock-in angle with field. Our results show that a variety of 
complex staircases are formed depending on the field orientation and strongly 
suggest that, at high temperatures, correlated structures dominate vortex 
pinning over random disorder in the whole angular range. 

Work partially supported by ANPCyT, Argentina, PICT 97 No.01120. 
A.S. and G.N. are members of CONICET. We acknowledge useful discus- 
sions with S. Grigera, F. de la Cruz, E. Osquiguil and D. Niebieskikwiat. 

Figure 1: Irreversible magnetization Mj as a function of the applied field angle 
for several fields H, at temperatures (a) 707^ and (b) 60K. Inset: blow 
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up of the H = 2T data near the c-axis for T = 60K (the units are the same 
of those in the main figure. 

Figure 2: Angle Qmax of the maximum in Mj(0) as a function of H for three 
temperatures. The sohd hues are fits to Eq.(l) (see text). Bars mark the 
width of the plateau. Inset: Mj(0) in the region of the plateau. 

Figure 3: Irreversible magnetization Mj versus field angle for three fields 
(curves are vertically displaced for clarity). Open symbols: data for > 0. 
Solid symbols: data for < 0, reflected with respect to the c-axis. Arrows 
indicate the angle Qgym beyond which the behavior is symmetric. The pro- 
cedure of reflection is sketched in the inset. 

Figure 4:Temperature dependence of Qgym (see fig. 3). The solid line is a 
guide to the eye. The sketches show possible vortex staircases for > Q^,. 



References 

[1] L. Civale et al, Phys. Rev. Lett. 67, 648 (1991). 

[2] W. K. Kwok et ai, Phys. Rev. Lett. 69, 3370 (1992) 

[3] I. V. Grigorieva, L. A. Gurevich, and L. Y. Vinninkov, Physica C 195, 
327 (1992). 

[4] W. K. Kwok et al, Phys. Rev. Lett. 64, 966 (1990); S. Fleshier et ai, 
Phys. Rev. B 47, 14448 (1993). 

[5] Y. H. Li, Phys. Rev. B 48, 6612 (1993). 

[6] L. Klein et al, Phys. Rev. B 47, 12349 (1993); B. Holzapfel et al., Phys. 
Rev. B 48, 600 (1993). 

[7] V. Hardy et al, Phys. Rev. Lett. 54, 656 (1996). 

[8] M. Oussena et ai, Phys. Rev. Lett. 76, 2559 (1996). 

[9] A. A. Zhukov et a/., Phys. Rev. B 56, 3481 (1997). 



[10] A. Herbsommer et al, Physica C 304, 112 (1998). 

[11] D. R. Nelson and V. M. Vinokur, Phys. Rev. Lett. 68, 2398 (1992); 
Phys.Rev.B. 48, 13060 (1993). 

[12] G. Blatter et al, Rev. Mod. Phys. 66, 1125 (1994). 

[13] F. de la Cruz, D. Lopez and G. Nieva et al., Pliilos. Mag. B 70, 773 
(1994). 

[14] S. Candia and L. Givale, submitted to Supercond. Sci. Technol.; D. M. 
Casa et al, unpublished. 

[15] A. A. Zhukov, G. K. Perkins, Y. V. Bugoslavsky, and A. D. Caplin, 
Phys. Rev. B 56, 2809 (1997). 

[16] L. Krusin-Elbaum, L. Givale, J. R. Thompson, and G. Feild, Phys. Rev. 
B 53, 11744 (1996). 

[17] G. Blatter, V. B. Geshkenbein, and A. I. Larkin, Phys. Rev. B 68, 875 
(1992). 

[18] E. H. Morre et al, Phys. Lett. A 233, 130 (1997). 



